function SAT_SST_Step_05_calculate_anomalies_for_scaled_SAT(P)

    % *********************************************************************
    % Load Scaled Air tempertures
    % *********************************************************************
    file_load = [SAT_SST_func_IO('data',P),'/Raw_SAT_SST_',P.date,'.mat'];
    load(file_load,'lon','lat');

    file_load = [SAT_SST_func_IO('data',P),'/Scaled_SAT_full_',P.date,'.mat'];
    load(file_load,'Scl_AT','Scl_AT_cowtan');

    % *********************************************************************
    % Call external function to connect stations and calculate anomalies
    % *********************************************************************
    [Scl_AT_anm,~] = SAT_SST_func_connect_stations(Scl_AT,lon,lat,...
                            P.do_gradual_match,P.distance,P.cal_anomalies);

    [Scl_AT_anm_cowtan,~] = SAT_SST_func_connect_stations(Scl_AT_cowtan,lon,lat,...
                            P.do_gradual_match,P.distance,P.cal_anomalies);

    % *********************************************************************
    % Save data
    % *********************************************************************
    file_save  = [SAT_SST_func_IO('data',P),'/Scaled_SAT_anm_',P.date,'.mat'];
    save(file_save,'Scl_AT_anm','Scl_AT_anm_cowtan','-v7.3');

end
